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Abstract. Consider a time-varying collection of n points on the positive real 
axis, modeled as Exponentials of n Browniaii motions whose drift vector at 
every time point is determined by the relative ranks of the coordinate processes 
at that time. If at each time point we divide the points by their sum, under 
suitable assumptions the rescaled point process converges to a stationary dis- 
tribution (depending on n and the vector of drifts) as time goes to infinity. 
This stationary distribution can be exactly computed using a recent result of 
Pal and Pitman. The model and the rescaled point process arc both central 
objects of study in models of equity markets introduced by Banner, Fcrnliolz, 
and Karatzas. In this paper, we look at the behavior of this point process 
under the stationary measure as n tends to infinity. Under a certain 'conti- 
nuity at the edge' condition on the drifts, we show that one of the following 
must happen: either (i) all points converge to 0, or (ii) the maximum goes to 
1 and the rest go to 0, or (iii) the processes converge in law to a non-trivial 
Poisson-Dirichlet distribution. The underlying idea of the proof is inspired by 
Talagrand's analysis of the low temperature phase of Derrida's Random En- 
ergy Model of spin glasses. The main result establishes a universality property 
for the BFK iikhIcIk and aids in explicit asymptotic computations using known 
results about the Poisson-Dirichlet law. 



1. Introduction 

Let 61,62, ■■■ ,Sn be n real constants, and consider the following system of sto- 
chastic differential equations: 

n 

(1) dXi{t) = ^ (5^1 {Xi{t) = dt + dWi{t), i = 1, 2, . . . , n. 

Here X(^^{t) > X(2)(t) > • . • > X(„)(t) are the coordinates of the process in the de- 
creasing order, and W = (VFi, W2, . . . , Wn) is an n-dimensional Brownian motion. 
The SDE models the movement of n particles as interax;ting Brownian motions 
such that at every time point, if we order the positions of the particles, then the ith 
ranked particle from the top gets a drift 5, for i = 1, 2, ... , n. As time evolves, the 
Brownian motions switch ranks and drifts, and hence their motion is determined by 
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such time dependent interactions. Such processes (or close relatives) have been con- 
sidered recently in several articles. We consider three major veins in the literature 
where these processes make their appearance. 

1. Stochastic Portfolio Theory. A detailed study of the solutions of SDE ([l]) 
in general was taken up in a paper by Banner, Fernholz, and Karatzas (BFK) [6]. 
These authors actually consider a more general class of SDEs than ([T|) in which the 
drifts as well as the volatilities of the different Brownian motions depend on their 
ranks. Using a delicate and elegant analysis of the local time of intersections of the 
different Brownian motions, the authors were able to demonstrate various ergodic 
properties of such processes. An important object in their study is the vector of 
spacings in the ordered process: 

(2) :=X(,)(i)-X(,+i)(t), l<i<n-l. 

Under suitable sufficient conditions the authors proved that the marginal distri- 
butions of the spacings converge to Exponential laws and left the question of joint 
convergence open. The joint law was identified and the BFK conditions were shown 
to be necessary and sufficient in Pal and Pitman [25 (see Theorem [l] below for the 
statement of the Pal-Pitman result). 

2. Burgers', McKean-Vlasov, and the granular media equations. Sznit- 
man in his articles on the propagation of chaos, [34 and [35 , considers n Brow- 
nian motions arranged in decreasing order. He proves that the ordered process 
can be thought of as a multidimensional Brownian motion reflected in the orthant 
{xi > X2 ■ ■ ■ > Xn} and discusses its relation with the solution of the Burgers' equa- 
tion. For more details on reflected Brownian motions in an orthant, see Varadhan 
and Williams ^37] and Williams [3^ . For more on stochastic particle system meth- 
ods in analyzing Burgers' equation, see Caldcroni and Pulvirenti [TT]. Note that 
the drifts in SDE ^ at any time point are a function of the empirical cumulative 
distribution function induced by the n Brownian particles on the line. In stochastic 
models related to McKean-Vlasov and the granular media PDE, the empirical CDF, 
for large n, is an approximation to a single probability measure (depending on t). 
Thus, ([l]) and similar systems often appear as approximations to SDEs whose drifts 
and volatilities depend on the marginal distribution of the particles. For example, 
in the context of McKean-Vlasov type equations, see Bossy and Talay [TU] and the 
recent article by BoUey, Guillin, and Villani [9]. For the granular media equation, 
see the article by Cattiaux, Guillin, and Malrieu [T2] which also lucidly describes 
the physics behind the modeling. Also see Jourdain [53] and Malrieu [57] for an 
integrated approach to the McKean-Vlasov and granular media differential equa- 
tions and related interacting diffusions. From similar perspectives Jourdain and 
Malrieu [25] considers SDE ([TJ with an increasing sequence of Si's and establishes 
joint convergence of the spacing system ^ to independent Exponentials as t goes 
to infinity. 

3. Interacting particle systems and queueing networks. The ordered par- 
ticle system of Brownian motions had been studied even earlier by Harris in |15j . 
Harris considers a countable collection of ordered Brownian motions with no drifts 
(i.e. 5i = 0), and analyzes invariant distributions for the spacings between the 
ordered processes. Arratia in [5] (see also [4]) observes that instead of Brownian 
motion, if one considers the exclusion process associated with the nearest neighbor 
random walk on Z, then the corresponding stationary system of spacings between 
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particles can be interpreted as a finite or infinite series of queues. We direct the 
reader to the seminal articles by Harrison and his coauthors for a background on 
systems of Brownian queues: |T7], [TS], [S], and [217. Baryshnikov ^7^ estab- 

lishes the connections between Brownian queues and GUE random matrices. For 
more references in this direction, see Pal and Pitman [28 . Also, see the related 
discrete time models in statistical physics studied by Ruzmaikina and Aizenman 
|32j . and Arguin and Aizenmann [2\. In particular see the article by Arguin |3] 
where he connects the Poisson-Dirichlet's with the competing particle models. 

The article by Pal and Pitman [28] considers solutions of SDE ([T|) as their cen- 
tral object of study. The number n can either be finite or countably infinite. They 
show the existence of a weak solution for SDE (jlj and the uniqueness in law (see 
Lemma 3 in [H]). One of the major results in the paper is a necessary and suffi- 
cient condition on the sequence {6i, . . . , (5„} such that the law of the spacing vector 
(Yi(i), Y2{t), . . . , Yn-i{t)) converges in total variation to a unique stationary distri- 
bution as time goes to infinity. We would like to caution that 28] uses the order 
statistic notation in which Xji) < X(^2) < • ■ • while describing SDE (jlJ, which is 
the reverse of the notation we found to be the most suitable for this paper. We 
state part of the result below following our notation in ([T]). 

Theorem 1 (Pal and Pitman [28 , Theorem 4). Consider the SDE in ([T]). For 

I < k < n, let 

k 

(3) ak 5^(4 - S,) 

i=l 

where 6^ is the average drift (Si + 62 + ■ ■ ■ + 5n)/n. 

For each fixed initial distribution of the n particle system with drifts 5i, the 
collection of laws of X(^i'^(t) — X(„-)(t) for t > is tight if and only if 

(4) ak > for all I < k < n - 1, 

in which case the following result holds: 
The distribution of the spacings system 

(r,(t) = X(,)(t) - X(,+i)(t), 1 < J < n - 1) 

at time t converges in total variation norm as t 00 to a unique stationary 
distribution which is that of independent Exponential variables Yj with rates 2a j, 
1 j n — 1. Moreover, the spacings system is reversible at equilibrium. 

Remark. The necessity of condition ^ is hinted by the following localization 
argument: for every k, one divides the collection of n particles as the top k particles 
and the bottom n~ k particles. If the entire collection has to remain together, the 
average drift of the top k particles must be less than than the average drift of the 
bottom n — k. A simple computation then implies (|4]). 

Now suppose we have a triangular array of drifts {Si (n) , 1 < i < n, n e N) , such 
that for every n the vector {5i{n) , S2{n) , . . . , Sn{n)) satisfies condition Note 
that the market weights (|9| are purely a function of the spacings, and thus have 
an induced stationary measure for every n. Thinking of the market weights under 
this stationary measure as a random point process in [0, 1], we are interested in the 
limit law of this point process as we take n to infinity. It turns out that any non- 
trivial limit law is a Poisson-Dirichlet distribution under 'universal conditions', and 
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exhibits phase transitions. The results are summarized in the following theorem. 
For the definition of the Poisson-Dirichlet family of distributions, see Section |3] 

Theorem 2. Suppose for every n (£N we have a sequence of constants (Si{n), i — 
1, 2, . . . , n), and consider n interacting particles satisfying the following SDE: 

n 

dX,{t) ^^Sj{n)l {X^{t) ^ dt + dW^{t), i^l,2,...,n. 

Here ^(i) > ^(2) 2i • • ■ 2i -^(n) f^f^ the ordered processes, and W is an n- 
dimensional Brownian motion. 

Suppose for every n e N, the drift sequence {Si{n), I < i < n) satisfies condi- 
tion Q. Then by Theorem^ the decreasing point process 

/ e-^u)(t) 

converges in distribution as t ^ oo. Let . . . , /i(„)(n)) denote a point process 

(of the market weights) drawn from the stationary law. Now let 5{n) — {5i{n) + 
... + 5n{n))/n and assume that there exists rj G [0,oo] such that for every fixed 
i > 1, 

(6) lim (S{n) — Si{n)) ~ rj. 

n — >oo 

Additionally, assume that 

(7) limsup max ((5(n) — (^^(n)) < rj. 

Then, as n tends to infinity, 

(i) if rj €z (0, 1/2), the point process generated by the market weights converges 
weakly to a Poisson-Dirichlet process with parameter 277. 

(ii) If r] €1 [1/2, cxd], all the market weights converge to zero in probability. 

(iii) 7/77 = 0, the market weights converge in probability to the sequence (1, 0,0,.. .), 
i.e., the largest weight goes to one in probability while the rest go to zero. 

In particular, converges weakly to a non-trivial distribution if and only if 

7? e (0,1/2). 

The conditions ^ and Q admit natural interpretations. Condition ^ is the 
'continuity at the edge' condition mentioned in the abstract. It means that the 
top few particles get almost the same drift. This is a weaker version of a general 
continuity condition that would demand that whenever two particles have almost 
the same rank, they would get almost the same drifts. The condition ^ can be 
understood as saying that the highest ranked particles get the minimum drifts, in a 
limiting sense. Although the conditions seem to be reasonable, it is still surprising 
that the conditions are also sharp, in the sense that Poisson-Dirichlet convergence 
may fail if either ([6]) or ([T]) fails to hold as shown by the counterexamples in Section 
|5] Finally we would like to mention that the two degenerate limits can be seen as 
limits of the Poisson-Dirichlet laws as the parameter 2ri tends to zero or one. 

Now, Theorem [2] says that /ii(n) ^ in probability if 77 > 1/2. It is natural to 
wonder about the rate of convergence. It turns out that some exact information 
may be extracted after putting in an additional 'Lipschitz' assumption on the drifts, 
even at the 'critical value' 77 = 1/2. 
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Theorem 3. Assume the same setup as in Theorem^ Let 

•q = lim ((5(n) — (5i(n)), 

n — *oo 

and instead of (|6| and assume that there exists constants C > and < 7 < 1 
such that whenever i < jn, we have 

(8) |5.(n)-(5i(7i)| < 

Then, if rj > 1/2, we have 

logfXi{n) 1 

y 1 m probability as n ~+ 00. 

log n 2rj 

If "q = 1/2, and we have the additional condition that 5{n) — 5i{n) = r] + 0{l/logn), 
then 

logMi(n) . 

> — 1 m probability as n 00. 

log log n 

Note that Theorem |2] describes a different fcind of hmit than the other possibihty 
of starting with a countably infinite collection of rank-dependent Brownian motions 
and describing their stationary behavior. To see the difference compare with the 
analysis of the Infinite Atlas model done in . 

The paper is organized as follows. Section [2] discusses the implications of the 
convergence to Poisson-Dirichlet law for the Banner-Fernholz-Karatzas models of 
equity markets. In particular, we compute the asymptotic laws of certain func- 
tional of the market weights by resorting to similar computations for the Poisson- 
Dirichlet law. They include the expected market entropy S = — J2i l^i log l^i ^^'^ the 
expected pth moment Dp = fif. Both these functions are considered in stochas- 
tic portfolio theory as measures of market diversity (see [13, page 31]). In Section |3] 
we introduce the Poisson-Dirichlet law and discuss its relation with certain Poisson 
processes on the real line. Section |4] contains proofs of all the major theorems in 
this paper. Finally, in Section |5] we discuss counterexamples demonstrating the 
tightness of the conditions in our theorems. 

2. Application in BFK and related models 

Fernholz in his 2002 book [T^ introduces solutions of ([T|) to model the time 
dynamics of the logarithm of the market capitalizations of different companies in 
an equity market. In other words, he considers a stock market with n companies 
whose total worth in stocks are coordinate-wise Exponentials of the solution of 
equation ([l]). A major objective of his work is to explain the following curious 
empirical fact regarding what are known as the market weights. The market weight 
of a company is defined to be the ratio of its market capitalization to the total 
market capitalization of all the companies, that is, the proportion of total market 
capital that belongs to a particular company. The market weight is a measure of 
the influence that the company exerts on the entire market and have been studied 
extensively in the economics literature over the decades for their peculiar structure 
and very stable behavior in time. For example see articles by Hashemi |21j . Ijiri 
and Simon Jovanovic and Simon and Bonini The structure in the 
market weights arises if one arranges the market weights data in decreasing order 
and plots their logarithms against the logarithm of their ranks. This log- log plot is 
referred to as the capital distribution curve. For example, if we consider Fernholz's 
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Figure 1. Capital distribution curves: 1929-1999 

model, and denote hy X a solution of SDE ([T|), the ordered market weights are 
given by the transformation 

The capital distribution curve will then be a plot of log^i versus logi. Empirically 
this plot exhibits nearly polynomial decay with decreasing ranks. See Figure 1 below 
(reproduced from |13p which shows capital distribution curves between 1929 and 
1999 for all the major US stock markets (NYSE, AMEX, and NASDAQ) combined. 
To explain this fact, Fernholz considers a special case of ([!]) called the Atlas model 
where (5„ > 0, and Si — for all i ^ n. He conjectured that the capital distribution 
curve obtained from the diffusion satisfying the SDE for the Atlas model converges 
to a stationary distribution under which it is roughly a straight line. 

It is interesting to note that the phase transition phenomenon described in The- 
orem [2] predicts some remarkable outcomes about models of stock market capital- 
izations as given in Banner, Fernholz, and Karatzas [6], [13], and [14] . 

As pointed out in [6], it is empirically observed that the companies with the 
few largest market weights get the least drift, and most movement seems to come 
from the bottom of the pile. This validates both our assumptions about the drift 
sequences made in ([6| and (|7]). Hence our results show the following universality 
property of the BFK models. When the number of companies in the market is 
large, one of the following must happen, either, the market share of each company 
is negligible with respect to the entire market, or the emergence of a 'superpower' 
which dominates the market, or a delicate equilibrium under which the market 
weights are arranged as a particular realization of a Poisson-Dirichlet process. These 
conclusions hold, no matter what specifically the values of the drifts are, as long as 
conditions ([6| and ([t]) are satisfied. The atoms of Poisson-Dirichlet are well-known 
to decay polynomially (the slope of the log-log plot for PD(q:) being roughly — 1/a) 
which validates Fernholz's conjecture. 
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The Atlas model (^) is an example where such a phase transition takes place. 

Example. Consider the Atlas model with n particles ([^, [H]): for i — 1,2, ... ,n, 
let 

dX,{t) = rjnl {X,{t) = X(„)(i)} dt + dW,{t), 77„ > 0, 

where W is an n-dimensional Brownian motion. The drift sequence is clearly given 
by Sn{n) = ?7„ and Si{n) = for i — 1,2, . . . ,n — 1. Thus, 6{n) — rin/n, and the 
drift sequence satisfies condition Q: 

ak = ^ {5{n) - 6,{n)) = ^ > 0, fc = 1, 2, . . . , n - 1. 

1=1 

Thus, by Theorem [T| there is a unique stationary distribution under which the 
spacings Y(^i^(t) = are independent Exponential random variables 

with rates 2ai. Also, note that if = lim„ rjn/n exists in [0, oo], then 

limsupmax (^(n) — (5i(n)) = 7], lim (S{n) ~ Si{n)'^ — r], i — 1,2,... 

Let be the stationary law of the spacing system in an n-particle Atlas model. 
It follows from Theorem [2] that as n tends to infinity, the point process of market 
weights under Qn either converges to a non-degenrate PD distribution, or to one 
of the two degenerate limit depending on the value of 77. Since the atoms of the 
Poisson-Dirichlet laws have polynomial decay, the Atlas model fits into the observed 
data in Figure 1. 

Moreover, the drift sequence trivially satisfies the Lipschitz condition ([s]) in The- 
orem |3] for 7 = 1/2. Hence, when 77 > 1/2 the market weights go to zero at a 
polynomial rate. Finally, suppose 

S{n) - = f]n/n = ^ + 0(l/logn). 

Then, the market weights go to zero at a logarithmic rate. In particular, this 
happens if 77,1 = 7i/2. 

We introduce another model, the one-dimensional gravity model, which satisfies 
such a transition. 

Example. For a parameter sequence r/„ > 0, let 

(10) dXiit) = '^y" sign (X J t) - X^(t))dt + dWi(t), i = l,2,...,n. 

n ^-^ 

Here {W\, W2, . . . , Wn) is an 77-dimensional Brownian motion. It is straightforward 
to verify that SDE ( 10 1 is a particular case of Q where the drift that the ith 
ranked particle gets is 5i — r]n(2i — n — l)/n. The more striking property of 
this model is the mutual attraction between particles which is due to the one- 
dimensional gravitational force kernel sign(-). One can think of the parameter 7y„ 
as the strength of gravitational pull between the particles. Note that the average 
drift for the gravity model is zero by symmetry. Moreover, the drift sequence 
satisfies condition (|4]): 



ak := Y{S{n) - Si{n)) = 



/fc(fc+l) fc(7l+ 1)^ _ 77„fc(77- fc) 

\ n n 
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which is positive for all k ^ l,2,...,n — 1. Thus, from Theorem [T] it follows 
that there is a unique stationary measure for the spacings (Yi(t), . . . , un- 
der which the coordinates are independent and the ith spacing is distributed as 
Exponential with rate 2ai. 

Now, suppose that rj :— lim„ ?]„ exists in [0,oo]. Then note that 

lim (S{n) - Si{n)) =77, i = 1, 2, . . . 

n — *oo 

lim max (S(n) — Si(n)) = lim77„ = rj. 

n l<i<n ' n 

Thus, the gravity model satisfies conditions ([6]), ([t]) for Theorem |2] Additionally, 
it follows that 

\5,{n) - 5^{n)\ < ^^'-^\ V n e N, 
n 

for some constant C depending on rj. Thus, the gravity model also satisfies condition 
([s]) for Theoremjs] We can thus conclude that all the conclusions of phase transition 
in Theorems [2] and [3] hold for the gravity model as 77 varies from zero to infinity. 

2.1. Some benefits of the Poisson-Dirichlet law. The limiting Poisson-Dirichlet 
law helps us to compute asymptotic laws of functions of the market weights under 
stationary distribution. Various distributional properties of PD(a), < a < 1, 
have been studied extensively, and we only use the most important ones for our 
purpose here. Many more identities together with their various applications can be 
found in [29] . 

The following theorem, whose proof can be found in Section |4j is a typical 
example of how such properties can be put to use in the context of BFK models. 

Theorem 4. For every n G N, consider the solution of the SDE ([ij with a drift 
sequence {5i{n)^52{n), . . . ,(5„(ri)}. Assume that the array of drifts satisfies condi- 
tions Q of Theorem^ and ^ and ^ of Theorem^ with rj G (0, 1/2). We have 
the following conclusions. 

(i) The asymptotic moments of the maximum market weight fJ-i{n) under the sta- 
tionary law is given by 



(11) Jim E(/.?(n))= tP-'e-'ij^^\t)dt, p > 0, 

where i^2ri{t) = 1 + 277 /o^(l ~ e^^^)x^'^'''~^dx. In particular, this determines the 
asymptotic law of the highest market weight. 

(ii) For p > 2ri, let Dp{n) = X]r=i l^ii''^)- Then the asymptotic expectation of Dp{n) 
is given by 

(in) Define the market entropy S{n) = ~ X^T^i Mi(^) ■ Then 

00 ^ 



RANK DEPENDENT BROWNIAN MOTIONS 



9 



3. Point processes and the Poisson-Dirichlet law 

In this section wc introduce the Poisson-Dirichlet law of random point processes 
and discuss the notion of weak convergence on the underlying space. 

Recall that a point process can be viewed as a random counting measure. For 
example, the Poisson point process with a cr-finite intensity measure fi, PPP(//), is 
defined to be a random point measure with the following properties. For any collec- 
tion of disjoint Borel subsets {Ei,E2, . . . ,Ek} oi the real line satisfying fJ,{Ei) < oo 
for all i, let N{Ei) be the number of elements of the PPP(/^) that are contained 
in Ei. Then the coordinates of the random vector {N{Ei), N{E2), ■ ■ ■ , N{Ek)) are 
independent and the ith coordinate is distributed as Poisson(/i(i?i)). If the ran- 
dom point process has a maximum atom almost surely, it is customary to order 
the atoms as a random decreasing sequence of real numbers. This happens, for 
example, when the intensity measure fi for PPP(/i) is such that oo) < oo. 

To discuss weak convergence of such random sequences it is necessary to define 
an underlying Polish space. The separable complete metric space (5, d) we use for 
this article is the space of all real decreasing sequences with the metric of pointwise 
convergence, i.e., 

(14) S = {x={x,), x,>X2>...}, d(x,y) = f;fc|^. 

i=l 

The fact that this is a Polish space and other properties can be found in page 45 
of |30j . From now on, our point processes will be random sequences taking values 
in the state space (5, d). Any finite sequences (xi, X2, ■ ■ ■ , Xn) will correspond to 
an infinite sequence {xi,X2, . . . , a;„, — oo, — oo, — oo, . . .) in (5, d). 

Here is a typical example of a point process which will be of use to us later 
on. Consider a sequence of independent standard Exponential random variables 
r]i, i — 1,2,.... Let »7i(n) > T]2{n) > ... > rjn^n) be the decreasingly ordered 
rearrangement among the first n values of the sequence. That is 

r]i{n) = max 77i, i]n{n) = min rji. 

l<'i<n l<i<n 

Lemma 5. The random point process s'^ = (?7i(n)— log n), i= 1,2, ... ,n, converges 
weakly to a Poisson point process on the line with intensity dfi* = e^^dx. 

Proof. Convergence of the extreme order statistics of iid data to a Poisson point 
process is a standard exercise in extreme value statistics. For example, see the 
book by Resnick |3Q| . In general, it follows from the fact that suppose, for every 
n > I, Xi, X2, . . . , Xn are n iid data with law P„. Further assume that there 
is a measure /i such that lim„_+oo "-^'n(^) = m(^) for all Borel subsets A of M. 
Then the random point measure induced by the n data points converge weakly 
to a Poisson Random measure with intensity /i. When the iid data comes from 
Exponential 1) distribution shifted by log 71, fJ.{A) = J^e^^dx — n*{A), and this 
proves the Lemma. □ 

The other point process important in this paper is a one parameter family of 
laws of point processes known as the Poisson-Dirichlet with parameter a, PD(a), 
for < a < 1. The PD family is fundamental and comes up in several contexts. It 
can also be seen as a part of a wider two parameter family of laws. For an excellent 
survey of old and new results we suggest the article by Pitman and Yor [25] . 
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PD(q;) can be defined in terms of a PPP with intensity measure diia = x^'^^^dx 
on (0,00). Let {Xi) be the decreasingly arranged atoms of PPP(/ic(). Then, by 
Proposition 10 in the point process 

^ — v-\oo Tr~ J * — 1 J 2, ... , 

is distributed as PT){a). We take this to be the definition of the PD family. It can 
be shown that each Vi is strictly between zero and one, and Vi = I. For more 
details, see [25]. 

Now the Poisson point process in Lemma [5] has a very fundamental connec- 
tion with the Poisson-Dirichlet (PD) process which is made clear through the next 
lemma. 

Lemma 6. Let {Yi} be a realization of a Poisson point process in decreasing order 
with intensity fj,* as in the last lemma. Then for any (3 > 1, the law of the point 
process generated by the sequence 

_ exp/gy, • _ 1 9 
is the Poisson-Dirichlet distribution PD{\/[3). 

Proof. From Proposition 5.2 in |301 page 121], it follows that the point process 
exp l3Yi is a Poisson point process on the real line of intensity [3~^x~^/^~^dx. Since 
/3 > 1, the result follows from the definition of PD(l//3) once we note that the 
constant in the intensity measure is a scaling parameter which gets cancelled 
when we take the ratio. □ 

The following theorem has been proved by Talagrand in [36' Lemma 1.2.1, page 
15] in the case of iid Gaussians. In Section |4] we follow more or less the same 
argument to prove it for Exponential random variables. The difference in the two 
arguments is merely due to the fact that the tails of Exponentials are fatter than 
those of Gaussians. 

Theorem 7. Consider again the ordered Exponential random variables (rjiin), i = 
1, 2, . . . , n) . For any /3 > I define the point process 

rr* _ ( exp{/3??,(n)} ■ _ 1 9 ^ 

Then, as n tends to infinity, U* converges weakly to the PD(\/[3) law. 

The difficulty in proving Theorem [7] lies in the fact that it does not follow from 
the convergence of the numerators of [/* to the Poisson point process of Lemma [s] 
Although Poisson-Dirichlet can be obtained from the limiting sequence as described 
by Lemma |6] the convergence of the point processes of the extreme statistics is in 
the sense of distributions which does not guarantee convergence of its sum or the 
ratio. The hard work lies in estimating the denominator so that one can show 
the convergence of the sequence of ratios directly to Poisson-Dirichlet. The same 
problem will be faced later, greatly amplified, when we prove Theorem |2] 

Weak convergence on the space (5, d) can sometimes be conveniently proved by 
suitable coupling thanks to the following version of Slutsky's theorem. 
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Lemma 8. Suppose on the same probability space we have two sequences of point 
processes X{n) = {Xi{n), i = 1,2, . . .} and X'{n) = i = 1, 2, . . .}, n G N. 

Suppose X(n) converges in law to X , and Xi{n) — X[{n) goes to zero in probability 
for every i, as n tends to infinity. Then X' {n) also converges weakly to X. 

Proof. The result follows from Slutsky's Theorem (Theorem 3.4 in [30j) once we 

show that d{X {n) X' (n)) 0. To achieve this, note that the random variables 
— Xj'(ri)| A 1 < 1 and converges to zero in probability by our hypothesis. A 
simple application of dominated convergence theorem imply that 



lim Ed(X(n),X'(n)) = lim eV 



\X.,{n)-X[{n)\M 



This shows that d{X {n) , X' {n)) converges to zero in probability and proves the 
lemma. □ 

Now, the support of the PD law is over sequences which have a special structure, 
namely, all the elements are non-negative and they add up to one. It will be easier 
for our arguments if we now restrict our sequence space to exactly such sequence 
with an apparently stronger metric. Let S' be the set of all sequences (xi) such 
that xi > a;2 > . . . > for all i and ^ 1- For any two sequences {xi) and {yi) 

we also consider the L^-norm between them 



(15) d'{x,y)^ 



El 



Weak convergence of measures on (5', d') is determined by functions which are 
continuous with respect to the L^-norm. In fact, we will, without loss of generality, 
consider functions which are Lipschitz with Lipschitz coefficient not more than one. 
That is to say, a sequence of probability measures /i„ on (S' , d') converges weakly 
to another probability measure /i on the same space if and only if 

(16) j fd^n J fdn, 

for all / which are Lipschitz with Lipschitz coefficients not more than one. 

Finally a note on the notation used in the following section. The Poisson-Dirchlet 
law is a probability law on (S',d'). For any test function /, we will denote its 
expectation with respect to any PD law by PD(/). The parameter of the PD law 
will be obvious from the context. 



4. Proofs 

4.1. Proof of Theorem [7} Let us recall the set-up: 771, 772, . . . is a sequence of iid 
Exponential one random variables, and 771 (n) > V^in-) ^ ■ • • > Vnin) be the ordered 
values of the first n random variables. For any /? > 1, we want to find the limiting 
law of the point process 

""lE"=iexp{/377,(n)}' ^ " 1' ^' • • ■ ' " I • 
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The first step is to fix a real number b, and define 

h,{n) := T],{n) - logn, u,{n) := e'^'*'^"^ 

(17) " " 

S{n) := ^Ui{n), Sb{n) := Ui(n)l{ft,i(n) > b}. 

i=l 1=1 

Let us denote Wi(n) = Uiin) / S{n), i = 1, 2, . . . , n, which are the coordinates of the 
point process J7*. We also define another point process depending on b: 



(18) u.,(6,n) 



0, if hi{n) < b, 

Ui{n) / Sb{n) , if hi{n) > b. 



The following lemma, exactly as in [351 Lemma 1.2.4, pg 17], estimates the difference 
between the two point processes. 

Lemma 9. 

S{n) - Sb{n) 



\wi{n) - Wi(6, n)| = 2- 



S{n) ■ 

Lemma 10. Given e > 0, we can find 6 G M such that 

(19) lim supP f 7 

Proof. Note that if S{n) < exp (3x for some a; G M, then rji — log n < x for all 
i — 1,2, ... ,n. The probability of this event is (1 — e~^/n) , which tends to 
exp(— exp(— a;)) as n tends to infinity. This proves that for any e, there is a 5 > 0, 
such that for all large enough n, one has P{S{n) < S) < e/2. 

We will show that given e, 5 as above, there is a 5 such that for all n large enough 

(20) P{S{n)-Sb{n)>eS)<^. 

It follows for that choice of b that for all n large enough we have 

S{n) 



In other words, (19 1 holds. 

To show ( 20 1 , note that by symmetry 



E(S'(n) - Sb{n)) = nE (e'^(''i-'°s«h{?7i - log?i < 6}) 



flog n+b 

= n^-'^ I e'^'-'e-'^dx 







[n''-iexp((/3-l)6)-l] 



eKP-i) 1 



/3-1 n/5-i(/3-l)' 

Since /3 > 1, we can choose 6 to be a large negative number such that for all n 
sufficiently large E(S'(n) — S'f,(n)) < €^5/2. Now (20 1 follows by Markov's inequafity. 

□ 
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Proof of Theorem^ The proof now follows exactly as in Talagrand's proof of The- 
orem 1.2.1 in page 19]. The basic idea being that from the last lemma it follows 
that the denominator in U* can be approximated by only about a Poisson many 
variables no matter what n is. The rest follows by pointwise convergence of these 
finitely many coordinates to the corresponding atoms of the Poisson point process 
which in turn approximates the Poisson-Dirichlet ratios. We skip the details. □ 

4.2. Proof of Theorem |2} The proof of Theorem |2] comes in several steps. The 
main idea is the following: compare the ordered process in stationarity with the 
point process of the order statistics of Exponential distribution as considered in 
Theorem [t] Condition ^ can be interpreted as saying that the top few atoms 
of both point processes have the same joint distribution up to a shift. The PD 
convergence requires exponentiation and taking a ratio as in Theorem [7] With 
the help of condition ^ we will show that when we exponentiate, the top few 
atoms gets so large compared to the rest that the denominator of the ratio in ([5| 
is basically a finite sum. This is the main technical difficulty in the proof which 
requires computations of suitable estimates. The PD convergence then follows from 
Theorem [T] 

The first step is the following lemma. 

Lemma 11. Suppose a triangular array of constants 1 < i < n—1), n>2, 

satisfies the following two conditions 



(21) limsup 



max 

l<i<n-l i 



< 1, and lim l3i{n) — i, i — 1,2, 3, 



Corresponding to this array of constants, let us take an array of Exponential ran- 
dom variables (Ci('^)j ^2('^); ■ • • j Crt-i('^)); n d N, all independent with E(^j(n)) = 
l/Pi{n). For each n G N and /3 > 0, let Un{0) he the point process with n coordi- 
nates given by 



t/„(/3)(j) = ^ 7 / ^, j = l,2,...,n-l, ano 



Unmn) 



1 



EL!exp(/3X;r=fc^C»(") 

1 

Er;!exp(/3x;r=fc 6(") 



Then 

(i) For any /3 > 1, the sequence Un{P) converges in distribution to the Poisson- 
Dirchlet law PD(1/ (3) as n oo. 

(ii) For any positive sequence {Sn} such that 5n < 1 and lim„^oo <5n = P ^ [Oj1]i 
the point process Un{6n) converges in law to unit mass on the sequence of 
all zeroes. In particular, J7„(/3) converges to the zero sequence when all the 
5n 's are identically equal to (3 for some /? < 1 . 

Let us first complete the proof of Theorem[2]using the above result. First suppose 
that 77 > in ([6]) and ([t]). In the notation used in Lemma 11 we define 



1 

(22) P,{n) = -Y,{5{n)~5,{n)), i^l,...,n-l. 
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Note that each Pi{n) > by Q. Thus, by assumptions (|6| and ([t]), we get 
(23 limsup max ^-^=limsup max -> < 1. 

n l<i<n-l i n I l<j<n-l i rj j 

Also, clearly 

lim 



= Z. 



Thus the array {Pi{n)} satisfies both the conditions in pTj ). Now, by Theorem[T] for 
the model with n particles, the joint law of the n spacings — X(^i^i-){t), i = 

1, 2, . . . , n— 1) under the stationary distribution is the same as that of a vector of in- 
dependent Exponentials (Yi(n), Y2{n), . . . , y„_i(n)) such that Yi(n) is Exponential 
with rate 

■i 

j=i 

which is the joint law of /2?y, . . . , i = 1, 2, . . . , n — 1). Thus, the decreasingly 

arranged market weights in the n particle system are in law given by 

exp(X(,)(n) - X(„)(n)) 



exp(X(j)(n) - X(„)(n)) 
' - J 7 — 12 n 

{n — 1 "1 ' — -L, z,, . . . , /I. 



From Lemma 11 it is immediate that if 77 G (0, 1/2), the point process of market 
weights converges weakly to the Poisson-Dirichlet law PD(277) as n tends to infinity, 
and if 77 > 1/2, the point process converges to zero. This proves conclusions (i) and 
(ii) in Theorem [2] 

Let us now prove part (iii) of Theorem [2j Assume that 77 = 0. For the n particle 
system, the largest market weight, fJ.i{n), is given by 

e-^(i)(") 1 



Eti e^«(") 1 + E:=^ exp{- ' 

where Yj(n) is Exponential with rate 2ai{n), where atiji) — ~ ^ji''^))- 

It is obvious that < 1. We will show that ^i(n)~^ ^ 1 in probability, 

which implies that ^ii{n) ^ 1 in probability too. Since all the market weights are 
nonnegative and add up to one, this forces the rest of them to converge to zero, 
and our theorem is proved. 

Now, recall that if F - Exp(A), then Ee^'^ is (1 + 1/A)"^ Thus 

n— 1 i 



l<E(l/Mi(n)) = l + ^nE(e-''^^"^) 
(24) 

^ ' n-l i ^ 

" "^^^n l + l/2a,(n)' 

1=1 7=1 ' ' 
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By assumptions ^ and Q, there exists no £ N such that 

max {S{n) — 6i{n)) < 2{5{n) — 5i{n)), V ?i > no. 

l<i<rL 

If we denote c„ = 5{n) — 5i{n), it follows from the definition of aj{n) that for all 
n > no we have aj{n) < 2jc„, and hence 

n— 1 i ^ n— 1 i ^ 

^ n 1 + i/2a,(n) - 51 n 1 + i/4jc„ ' 

Let 

'n-=i(l + l/4jc„)'\ ^<n, 
0, otherwise. 

Without loss of generality, we can assume, by (|6|, that c„ < 1/8 for all n > no. 
Thus 

^ n (1 + n = 2 (rfW - 

where C is some positive constant. Thus, qi{n) is dominated by a summable series 
for all n > no. Also note that, since c„ ^ by ([6]), we get lim„_>oo Qiin) = 0, V i. 
Hence we can apply the dominated convergence theorem to get 

oo 

limsup|E(l/^i(n)) - 1| < lim ^g,(n) = 0. 

n — >co ^ . T 

2 — 1 

Thus, we have shown that /ij^^ (n) converges to one in probability, and this completes 
the proof of Theorem [2] 

Now let us begin the proof of Lemma [TT] On a possibly enlarged sample space, 
consider an independent sequence {771, 772, . . .} of iid Exponential random variables 
with mean one. Let »7i(n) > rj2{n) > ... > r]nin) denote the first n random 
variables {rji, . . . ,r]n} arranged in decreasing order. As considered in Lemma [s] let 
£* denote the point process given by 

(25) £* = (77i(n) -logn, 7= 1,2, ...,n), n S N. 

For any 6 G N, define a new sequence of points by 



(26) l^'(n) 



?7j(n), j >b 

^7;b+i(n) + Et,6(n), 1<J<6. 

Let £„(&) be the new point process £„(&) = {Y-{n) — logn) , i = 1,2, We prove 

below that, for every 6 e N, as n tends to infinity, the law of en{b) converges to the 
same limiting law as that of £* as described in Lemma [Sj 

Lemma 12. For every fixed b £ N, as n tends to infinity, the random point process 
Sn(b) converges to a Poisson point process on the line with intensity dfj,* — e^^dx. 
Moreover, for every (3 > I, the normalized Exponential point process 

u (h\ / exp{/3y,(n)} 1 
[Ei=iexp{/3yj(n)} J 

converges in law to the Poisson- Dirichlet distribution with parameter 1//3. 
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Proof. Our main tool will be to apply Lemma |8] Let us recall the Renyi represe 
tation (131 page 114]): 



(28) 



(n) = r;(n) :=^?7,A, l<j<n, n G N, 



where — denotes equality in law. Also, trivially, the joint distribution of the se- 
quence {S,i{ri), 1 < i < n — 1} is the same as that of {rii/l3i{n), 1 < i < n — 1}. 



Thus, analogous to (26 1, it is clear that we can suitably redefine the sequence 
{iji} so that 



(29) 



Y;{n) 



- J 



For the rest of this argument we refer to e* as the point process given by {Y* (n) — 
logn, I < i < n). As before, e„ is the sequence of points {Y^{n) — logn, I < i < 
n). 

Recall that by Lemma |5] e* converges in law to a Poisson point process with 
intensity d^* = e"^c?x on the entire real line. To establish the same limit for the 
point process e„ we will verify the assumptions of Lemma [8j 

Now, by our construction, the difference between the jth coordinates of the two 
point processes |e*j(j) — £n(j)| is equal to 



\Y;{n)-Y;{n)\ < . , 



0, ifj>6, 



if 1 < J < fo, 



And thus 



E|l^*(n) - Y,'(n)| < 



b 

E 



1=3 



Now, the right hand side of the last equation goes to zero as n tends to infinity by 



the second condition in (21 1. It follows that the zth coordinate of e* — £„ goes to 
zero in probability, and the first part of the lemma follows by Lemma |8] 

Now, for any /3 > 1, by the continuous mapping theorem it follows that the 
difference between any coordinate of the point processes (exp(/3y/(n)), i < n) and 
{exp{(3Y* (n)) , i < n) goes to zero in probability. Also, since the two sequences 
differ only at finitely many coordinates, it follows that 



0. 



By another application of the continuous mapping theorem (see, e.g., |30| Theorem 
3.1, page 42]) we can conclude that the point process 



(30) 



exp{/3r/(n)} 



E;=iexp{/3i;'(n)}' 
converges to the same limiting law as the process 

exp{/3y;(n)} 



Y.UeMPYnn)y 



1,2,. 



1,2,. 
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By Theorem [7] and the equahty in law (28), we get that U*{b) converges to a 
PD(l//3) law as n tends to infinity. This completes the proof of the lemma. □ 



Define 
(31) 



1 < j < n - 1, 
i = n. 



Next, we will show that as h tends to infinity, the limiting point process for Un(b) 
is the same as the limiting point process (depending on the parameter /3) 



(32) 



Un 



exp/?X,(n) 
E?=iexp/3Xj(n)' 



1,2,. 



This will prove our objective that the limiting distribution for the point process 
Un is given by Poisson-Dirchlet with parameter 1//3. We start with the following 
lemma. 

Lemma 13. Let {rji, i ~ 1,2, . . .} be a sequence of independent Exponential one 
random variables. For (3 > I, let us define 



10, I = n. 

and let 

Let Un{b) be defined as in ( |30[ ). Then 

UmsupEd' (Un,Un{b)) 



I <i<n-l, 



E;.iexp/3Xj(n)^ 



z = 1,2,. 



< 



C 



where C is a constant depending on (3 and d' is the distance defined in (15 1. 

Proof. Note that, the sequence {X[{n), i = 1, 2, . . . , n} have the same joint law as 
the original sequence {Xi{n), i = 1, 2, . . . , n} in ( |31[ ). Thus, the new definition of 
the point process C/„ is consistent with the old one in ( [32| . 
Now, suppressing the dependence of b and /3, let 

exp(3Xl{n) ~ exp/3r/(n) 



C/.(n) 



E,exp/3X'(r 



E,exp/3y'(r 



Then, by definition, the difference J^^^i lUiin) — Ui{n)\ is equal to 



E 



exp/3(X,'(n)-X((n)) 



exp/3(y/(n)-y/(7i)) 



E, exp/3(X'(n) - X[{n)) exp/3(Y,'(n) - r/(n)) 



(33) 



i<b+l 

E 

i>b+l 



exp(3{Xi{n)-X[{n)) 



exp/3(r/(n)-r/(n)) 



exp {/3(Xj(n) - X[{n))} exp {P{Y^{n) ~ Y({n))} 

exp/3(XKn) - X[{n)) exp/3(y/(n) - r/(n)) 



exp/3(Xj(n) - X[{n)) exp/3(r;(n) - r/(n)) 
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Let Qn, Rm and Sn be defined (all depending on b and [3) as the partial sums 



i>b+l 



i>b+l 



Sn= Y exp{/3(X,'(n) - X^n))} = ^ cxp J -/3 ^ r;fe//3fe(n) I 

i<b+l i<b+l I fc<i J 



Y exp{/5(y/(n)-r/(n))}. 

i<6+l 



One can rewrite the first term on the RHS of ( 33 1 as 



i<b+l I k<i ) 



1 



1 



Sn 



Qn Sn Rn Sn 
Gn Rn 



Qn + Sn Rn + Sn 
Qn Rn 



]n + Sn Rn + Sn 



Qn + Sn Rn + "^n 

The final inequality follows by noting that Sn is greater than one. The rest of the 

exp/3(y/(n)-y/(n)) 



sum on the RHS of ( 33 1 is bounded above by 
exp/3(X;(n)~X(H) 

i> 



E 



E 



exp/3(Xj(n) - exp/3(r;(n) - y/(n)) 

Qn . Rn 



Sn Rn + Sn 



— Qn Rn ■ 



Hence the expected distance E d'{U{n), U{n)) is bounded above by 2E(Q„ + Rn). 
The proof is now completed by applying the next lemma. □ 



Lemma 14. 



limsupEQra < 



Ci 



limsupEi?„ < 



Ci 



where Ci is a constant depending on (3. 

Proof. We will first find a bound on the expected value of Qn . Note that 
Qn= Y ^^P"^ -/3E'?'=/^'=(") f- 

i>f)+l I /c<j J 

where {rjk, k = l,2,...n — 1} are a sequence of Exponential(l) and the array 



{Pk{n)^ fc = 1, 2, . . . , n — 1} satisfies condition (21 1 



EQ„ = EexpJ -/3^?7fc//3fe(n) Ie ^ exp J -/3 ^ r?fc//3fc(n) I 

[ k<b J i>6+l I 6<fe<z J 

1 " 1 

" n i + /3//3fc(„) E n l + /3//3fc(n) 
= Fn{h)-Gn(h) (say). 
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Now, since /? > 1, by (21 1, there exists uq £ N such that 

sup max Pi{n)/i < \f]i. 



n>no i^'i*^ 



In other words, for all n > no, and for all i = 1, 2, . . . , n — 1, we have < \f^i. 

Thus, for all n > no, we have 



4>6+l b<k<i I- I f I 

6<fc<i ^ i>h- 

(z-l)!r(6+ 



E 



By Stirling's approximation to the Gamma function, we can deduce (see, e.g., [TJ 
page 76, eq. 6.1.2]) that there exists a constant C depending on /3 such that 

<4=, V.>1. 



r(i + V?) - 

and thus one gets the bound 

limsup G„(o) < C > * 



^^r(6+y^+i) /-"^^ ^ ^,r(6 + v^ + i) 



where C" is another constant possibly depending on /?. 
Similarly, one can bound i^„(6) by 



lim^supF„(6) = nrTW 



< 



fe<b 

6! 



fe<,^ + V?/fc (V^+i)(V3+2)...(V^ + 6)- 



Thus, the product E(5„(6) can be bounded by 

„r(& + V^+i) 6! 



EQ„(6) < C 



< 



h\b^-^ (V?+i)(V^+2)...(v;g + 6) 

cT(v:g + i) Ci 



57?- 1 6v^-i' 

Here Ci is another constant depending on (3. This proves the first part of the 
lemma. 

The second part of bounding R(n) follows if we define the rates of the random 
variables i^/(n) as an array 



(3k{n) , l<k<b, 

k , b < k < n ~ 1. 



Then the array /3^.(n) also satisfy (21 1 and hence from the previous argument it 
follows that Ei?„(5) < Cib^^^ , and the lemma is proved. □ 
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We are now ready to finish the first part of Lemma [TT] We start with a test 
fimction on the space {S',d'), introduced in (15 1. That is to say a function / : 
5' — > M which is Lipschitz with Lipschitz coefficient one. Let PD{f) denote the 
expectation of / with respect to the PD(l//3) law. Now, by triangle inequality, 



(34) 



limsup |E/(C/„) - PD{f)\ < limsup E/([/„) - E/([/„) 

n — >oo n 

+ limsup Ef{Un)-PDif) 

n 

< limsupE d'(C/„, C7„) + limsup E/(C/„) - PD{f) 

n n 

c 



(/ — Lipschitz), 



< 



0. 



The first limit above is by Lemma [13] and the second limit is zero by Lemma |12[ 



We can now take h to infinity in (34 1, and since /? > 1, we conclude that 
lim„^oo E/([/„) = PD(f). Since this is true for all Lipschitz continuous func- 
tions, the conclusion is now evident from the standard theory of weak convergence. 



This completes the proof of part (i) of Lemma 11 



Let us now prove the second assertion in Lemma [TT] By Lemma [8] this is equiv- 
alent to proving that if < 5„ < 1, and lim„ Sn = f3 < 1, then each coordinate 
of the point process Un{Sn) converges in probability to zero. Obviously, it suffices 
to show that the first element of C/„((5„), which is also the maximum among the 
elements, goes to zero in probability. This will be shown in the rest of this section. 
Throughout the rest of this section, the parameter {(5„} will be a sequence such 
that each 6n is positive and less than or equal to one and lim„ Sn — f3 [0, 1]. 

As usual, let /ii(n.) denote the largest element in the random sequence C/„((5„). 
Then, we have 



(35) 



Mi(") = 



The following lemma, which proof wiU follow later, is our key tool in proving that 
^i(n) 0. 

Lemma 15. Suppose Vi, . . . ,Vk are independent Exponential random variables, 
with V, - Exp{Xi). Let 9^ = l/A^ = E{Vi). Let 

= 1 + e-Vi + e-'yi+v^) + . . . + e-(^i+-+^if) ' 
Let us also define the quantities 



(37) 



1 



K 
i=l 



We have 



(38) 
(39) 
(40) 



EM 

E (log fj.v - log p.vf 



> 
< 
< 



e fly, 



and 

Moreover, 
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Let us now complete the proof of Lemma [TT] with the aid of the last lemma. Fix 
some K > 2 and define 



(41) 



A*i K{n) :— ^— > ^i(n). 

1 + E,=2 exp{-(S„ ELi 



In the following, we will simply write /u„ and pLK,n instead of fJ-i{n) and /ii,if(n) 
respectively. For fixed n, if we take 9i(n) = Sn/Pi{n), and consider Vi as Exponential 
with mean 6i{n) for i — 1, . . . ,K, then it is clear that fiy has the same law fiK,n- 
Analogously, let 



1 + Ef=2exp{-ES^.H} 
Thus, for every e > 0, we can apply Lemma [T5| to get the Markov's bound 



(42) 



P(M„ > .) < P(MK,„ > .) < ^ < C^. 



Here C is the constant 4e^/*, and a (depends on K and n) is given by 



K 



Sn 



It follows from condition (21 1, and since (5„ < 1, that as n grows to infinity, while 



keeping K fixed, a has a limit given by 

aK — lim a^{K, n) = (3"^ 'S^ 



K 



i=l 



Thus, taking limit as n tends to infinity while keeping K fixed in (42 1, we obtain 

l/2aK 



(43) 



lim sup P(/i„ > e) < C ( -limsup/i/j- 



The RHS of (43 1 will be interpreted as zero if /3 is zero and limsup„ iJ,K,n < ^■ 
Suppose now that f3 > 0. In the next lemma. Lemma |16| we will show that 



(44) 



lim limsup /ij^ „ = 0. 

K >00 ji — 

^2^2 



Since ax grows to a finite limit /3 tt /6 as K tends to infinity, we take a further 
limit in ( 43 ) to get 



lim sup P(/i„ > e) < lim C -lim sup /iA' 



C - lim lim sup rifr , 
^ e K „ ' 



l/2aK 



0. 



This proves that /i„ goes to zero in probability. 

If /? = 0, for every e > 0, Lemma [16] shows that there exists K such that 

lim sup ]lK, n < e. 



We apply (43 1 to this K and obtain that limsup„ P(^„ > e) = 0. Hence we have 
established the second claim in Lemma [TT] completing its proof. 
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Lemma 16. 



lim limsup fiK,n — 0. 



n — ►oo 

Proof. We can write down fLx exactly in terms of the 9i{n) — Sn/Pi{n) as 

1 



As we take n to infinity keeping K fixed, it follows from condition (21 ) that 

n 

< 



< 



1 + 6-/5/1 g-(/3/l+/3/2) _\ ^ g-(/3/l+/3/2+-+/3//f) 

const. 

l + e-l^ + e-(/3iog2) ^ g-(/3iogK) 

const. 



Note that since /3 < 1 , the denominator of the last expression above goes to infinity 
as K tends to infinity. Thus 

1 1 



lim /iif = const, lim 1 + — ^ + . . . + — 77 = 0. 
This proves the lemma. □ 



Proof of Lemma 15 First, by Jensen's inequality, we have 

Now, by independence, S(e-(^i+ -+^')) = ni=i(l + ^'i)"^- Thus, 

(45) ^(a^^)^ (i + EITt^ 

Now, the function log(l + x) — {x ~ x"^) can be easily verified to be zero at zero, 
and increasing on the positive half-line. Thus, 

1 



It follows that 



> 



■ - = e Ml/. 

l + ELexp{-EU^J 



This proves (38 1. 



Next, consider the martingale 
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For the non- negative submartingale exp(Mi / Aa) and any p > 1, we apply Doob's 
i^'-inequality (|3I1 page 54]) to get 



(47) 



E ( max eP^^'/"^ ] < 

, l<i<K 



p-1 



PMk /4cr 



^ g-0./2o 



Now, by definition, 9i < a for every i. Thus, with p = 2, we get 

1=1 1=1 " 2a 

Now, as in the case of the upper bound, it is straightforward to verify that for any 
X € [0, 1/2], 1 — a; > e^^^^ . Thus, plugging in this inequality in the last expression 
we get 

(48) £;(e*^-/2-) < e^' "'/^'^^ = e'/\ 

Finally, note that for any two positive sequences {xi} and {yi} one has 

Xi+ X2 + ■ ■ ■ + Xii ^ Xi 

< max — . 

yi + y2 + ■■■ + yK i<«<^f Vi 

We use the above inequality to obtain 

fJ-V 



l + e""! H he^^ 

1 + e-^i + • • • + e-(^i+-+^'^) 



1 + 



(49) 



71 T _ 

< max e ' • uy. 

i<i<if 



Combining (|47|), (|48|), and (|49|), we get 

1/2(7 



l/4-l/2<T 



which proves ( 39 1 . 



Next, from M9k observe that 



(50) |£;(log^v) -logMvl < E( max \M,\) < 2[E{MI)Y/^ = 2a. 

Now recall that the Exponential distribution satisfies the Poincare inequality with 
Poincare constant 4 (see e.g. [8j, p. 2). That is, for any function f{Vi, . . . , Vk), 

(51) Var(/(l^i, ...,Vk))< 4E ^fi?(^(V;)' ' 



If we define 



f{xi,X2,.-.,XK) = - log 



K 

1 + ^ exp ■ 

i=i 



note that f{Vi,V2, ■ ■ ■ , Vk) — log^y. The partial derivatives of / are given by 



a/ 

5a;,; 



J2f=^ exp{-ELi2;fc} 
1 + Ef=i exp {- J2i=i } 



e[o,i]. 



Thus, a variance upper bound follows from (51 ): Var(log/iy) < 4(7^. Combining 
the last step with equation ( 50 1 , we get 

^(log/iy-log/iv)2<8a2. 
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This completes the proof of the theorem. 
4.3. Proof of Theorem [H 

Proof of Theorem^ In accordance to the set-up in Lemma 15 we define 

i 

X,{n)^2Y,iS{n)~5,{n)), 



□ 



and 



i{n) 



Aj(n) 



i — 1, . . . , n — 1. 



Now let Yi(n), . . . , y„_i(n) denote the successive increments in the stationary dis- 
tribution (as in Theorem |2|. Then Yi{n) has an Exponential distribution with rate 
Xi{n). Let 7' = {ri/2C) A 7, and define 

fJ-'iH ■■= r— , and 



i + El2'i"^exp(-E;=ii;H)' 

1 



i+Eiirexp(-E;.i^.H) 

Now let 

e„ = |2(5(n)-5i(n))-2?7|, 

and assume that n is large enough to ensure that e„ < rj/2. Since |(5i(ri) — Si{n)\ < 
C{i — ^)/n, it follows that 



^ i 1 

\K{n) - 2i'i]\ < eni + 2C V ^— < e„i + — , i = 1, 

n n 

If i < 7'n, it follows in particular that 

(52) X,{n) > 2iT] ~ e„i - Ci{-q/2C) > f]i. 

Thus, for n sufliciently large and i < "f'n, 



, n — 1. 



(53) 

Now let 



1 

2ir] 



Xi{n) — 2irj 
2i7]Xi{n) 



2ry 



< 



C 
2ri'^n 



Then for n sufficiently large and i < 7'n, 

' ' 1 a 

^ 0j(n) < a+ ^ ^ + < a+ logi + L(77), 

where L(ri) does not depend on n. Now, if r/ > 1 /2 then for all sufficiently large n, 
a+ < {l + l/2r))/2 < 1. Thus, 



(54) 



1 + e-^(") i-"" 



where k^(?7) is a constant free of n. Another application of (53 1 shows that 
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where a„ = (I/277) — e„ and k (rj) is another constant depending only on rj. Since 
hma^ = hma^ = 1/27], this shows that when rj > 1/2, 

hm l^MM^l.i. 
n^oo log n 2?7 



Again, by (52 1 it follows that 

(55) limsup ^ Oi{rLf < 00. 



We can now use Lemma |15I to conclude that 

log^'i(n) p 1 



1. 



n /^ 1 M 1 , SLl7'r»j+iexp(~Ej^i^j(")) ^^ 

I log Ml (n) - log Ml (n) I = log 1 + . ^, .. ) 

V 1 + >:,-i exp(- V_-, YAn)) J 



log n 2r] 

Finally, note that 

E7 =lYn\+i exp(~ E]=i Yjjn)) ' 

El2'i"Jexp(-E;.iiS-N). 

^ EmV«J+i exp(- Ej^i Yjjn)) 
- l + ElZ'i"^exp(-E;=ii^,H) 

^ Er= L7'nj +1 exp(- Ej-^i y?- ('^)) 

^ n— 1 y i 

I-' -I i=L7'riJ+l ^ j=L7'nJ+l 

^ n - - 1 ^ 1 - 7' 

~ [7'nJ +1 ~ 7' 

This completes the proof of the Theorem in the case 77 > 1/2. 

Now, in the case rj — 1/2, everything up to the first inequality in (54) is still 
valid. However, the second inequality does not hold since lima^ = 1. Instead, we 
use the inequality 

[7'nj Lt'jiJ 

> > Kn-'" logn 

i=l 1=1 

and the assumption that e„ — 0(1/ log rt) to conclude that 

fl[{n) < 



log n ' 

where K is some constant that does not depend on n. The inequality in the oppo- 
site direction follows similarly, possibly with a different constant. Inequality (55) 
continues to hold without any problem, and so does the subsequent argument. 
Combining, it follows that 

log Ml W _^ 
log log n 

This completes the proof. □ 
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4.4. Proof of Theorem IH 

Proof of Theorem^ The central argument in this proof is to apply continuous 
mapping theorem to the weak convergence result in Theorem [2] Let Qn be the 
stationary law of the market weights {/ii(ri), ^2("), ■ ■ ■ , fJ-nin)} arranged as a de- 
creasing sequence on the line. Then, by Theorem |2] Q„ converges weakly to the 
PD(277) law as n tends to infinity. Let {ct, i = 1, 2, . . .} denote a random sequence 
whose law is PD(2?7). 

To prove (i), note that the function x i— > is a continuous map in the space of 
sequences {S',d') (see (15l). It follows that Ati(n) converges weakly to the law of 
c^, and hence 



lim E^''(,^) = Ec? = — / tP-'e''ij^^\t)dt, 
^2^(t) = l + 2r] [ (1 - e-*")a;-2'?-idx, 



which proves (111. The expressions for the moment of the coordinates of a PD 
random sequence can be found in j29l Proposition 17] (for 9 — 0,a — 2ri,n = 1). 

Proving (ii) is similar. Again, by the continuous mapping theorem we know that 
for each i, 

lim E/if (n) = Ecf . 

n — *oo 

Now, for any bounded function / on [0, 1], the following identity that can be found 
in [291 page 858]: 

°° 1 r'^ (^ - 7;'|2'?-i 



If we take f{x) — x^, as before, it is easy to verify that 

r(p ~ 2r,) 

r(p)r(i - 27?) ■ 



(56) E^cf 



1=1 



Now, by the dominated convergence theorem, it is clear that lim„ E^ fJ-ii^^) equals 
the above expression if we can find a function ip : Z^. — > M such that (i) V'(*) < 

oo, and (ii) for all sufficiently large n, E/i^(n) < ip^i) for all i. 
Now, in the notation of Theorem [2l we have 

where lj(n)'s are independent Exponentials with rate 2a j{n) from Q. Now fix 
p' G (2-1], p). Let r„ = maxi<i<„((5(n) — Si{n)). Since lim„ r„ ~ 2rj, we can assume 
that n is sufficiently large to ensure that r„ < p' ■ Then 

i — 1 i—1 

E(uf (n)) < TT -, < TT -, =: iPii). 

;}i 1 + (") " 1 1 1 + pI^p 

It is not difficult to see that for large «, is comparable with I^p/p . Since p' < p, 
it follows that J2 V'(*) < and this completes the proof of part (ii). 
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For part (iii), let us begin with the definitions 

n oo 

(57) /i„(p) = E^/.f(n), /i(p) = E^cf, p>27j. 

i=l i=l 

By part (ii), hn{p) converges to /i(p) as n tends to infinity pointwise for every 
p > 2r]. Also, note that {hn} and h are differentiable convex functions in (277,00). 
It follows that (see, e.g., [221 Proposition 4.3.4]) their derivatives also converge 
pointwise in the relative interior (277,00). That is to say, lim„^oo ^Ii(p) — h'{p). 
Taking derivatives inside the expectation at 73 = 1 in (|57| , we get 



lim E V/i,(?^)log7ii(n) = h' {I) = E V q log i 



Now, to evaluate the last expression, we use the expression of h from (56): 

'V'{p-2ij) T{p-2r^)V'{p)- 



h'{p) = 



1 

r(i- 277) 



T^{p) 



Thus, 
(58) 



E^cJogCi = — 



r'(l-2r;) r'(i) 



(1 - 277) r(i)- 



It is known that (see, for example [H 6.3.2, 6.3.16, page 79]) 



ni) 

r(i) 



-7, 



r'(l-2ry) 

r(i - 27,) 



-7-2?7^ 



k=l 



k{k -2t])' 



where 7 is the Euler constant. Plugging in these expressions in (58 1, we get 

n ' 00 ^ 

^fii{n) log fi^{n) =2?7^^ 

i=l J 



lim E 

n — *oo 



fc=l 



k{k - 2r,) 



which proves the result. 



□ 



5. Counterexamples 

We now provide two examples where Theorem [2] fails due to the absence of one 
of the two conditions ^ and Q. 

Example. Let us consider SDE ([ij for n particles where 

(5i(7i) = -1/4, (5,(n)=0, z = 2,3,... 

This is like the Atlas model with parameter 1/4 except that the push comes from 
the top. Thus 

k 



(59) 5{n) 



1 

in 



> ak{n) ^^{5{n) - 5i{n)) 



k 
in 



n ~ k 
in ' 



which is positive for 1 < fc < n — 1. Thus, the drift sequence satisfy condition (|4]). 
Additionally the drift sequence satisfy ([t]) for 77 = 1/4, since 

lim max {5{n) - ^,(n)) = - - lim ^ = -. ^ (0, 1/2). 

However, the drift sequence does not satisfy ([t]). 
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It is easy to see that the market weights for this model do not converge to any 
Poisson-Dirichlet law. In fact, consider the following ratios 

/Lti(n) ' ^i2{n) 

where ^i(n) is Exponential with rate ai{n) for i = 1,2. From ( |59[ ), one can verify 
that lim„_+oo cn ('T^) = linin^oo Q^2(") = 1/4. Thus, it follows that the asymptotic 
laws of fi2{n) / fii(n) and ^^{n)/ ^2{n) are the same. However, if {Vi, V2, V3, . . .) is a 
random sequence from any PD distribution, it is known that (see Proposition 8 in 
[29] ) the laws of V2/V1 and V3/V2 are different Beta random variables. This negates 
the possibility that the market weights converge weakly to any PD law. 

Example. For the second example, consider any 77 e (0, 1/2). Let /3 — 4(1 — 77). 
For each n, define Si{n), i — 1, . . . ,n as follows: 

l<i< [n''J , 

Si{n) = { -/3, [n''J + l<i< [n/2\ , and 

i > [n/2\ . 

Then, by symmetry, S{n) = and the drifts satisfy the condition Q. 

It is also clear that the drifts satisfy (|6]), but ^ does not hold. We shall show 
that /ii(n) even though rj £ (0, 1/2). 

Let Yi (n) , i=l,...,n— Ibe the successive spacings under the stationary law 
as in Theorem [1] Then Yi{n) is Exponential with rate 

i 

K{n)^-2j2Sj{n). 

Note that by our construction, \n-i+i = ^i- Again, since /? > 77, it follows that for 
1 <i < 2[n''J we have 

(60) Xi{n) > 2rii. 

If 2\n^\ <i< Ln/2J , then i - [71^ \ >i/2, and we get 

(61) K{n) = 2r][n^\ +2(3{i- K'J) > (ii. 

Now let 6i{n) — l/\i{n). From the above observations and the summability of 
X^i^i it is clear that X)"^/ (") can be bounded by a constant that does not 
depend on n. So by Lemma |15[ we can conclude that ^i{n) — )■ in this model 
provided we can show that fii{n) — > 0, where 

^'^"^ ^ i + EL"/exp(-E;=i^,H)' 

Note that 

- / ^ 1 

/^l('^) < i — 7T\ ■ 

^''^ exp(E?ir;j^.(-)) 
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By (60), we get 



log n 
2 

for some constant C that does not depend on n. Again, for 2[n''J + 1 < i < [n/2\ , 
inequality (61 1 gives 



(63) E^.(")^^E 7^^ + ^ 



j=2[n')J+l j=2Ln'Jj+l-^ 

where C is again a constant that does not vary with n. Thus, 

L"/2J L"/2J 
(64) i=2LniJ+l j=2Ln''J+l 



where C" is some other constant. Combining the inequalities (62), (63), and (64), 
and using the relation = 4(1 — 77), we get 

This shows that /ii(n) ^ in probability. This shows that market weights do not 
converge weakly to any PD law inspite of condition ^ holding with 77 S (0, 1/2). 
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